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ABSTRACT 



We have considered non-conformal fluid dynamics whose gravity dual is a certain Einstein dilaton 
system with Liouville type dilaton potential, characterized by an intrinsic parameter n. We have 
discussed the Hawking-Page transition in this framework using hard-wall model and it turns out that 
the critical temperature of the Hawking-Page transition encapsulates a non-trivial dependence on 77. 
We also obtained transport coefficients such as AC conductivity, shear viscosity and diffusion constant 
in the hydrodynamic limit, which show non-trivial n dependent deviations from those in conformal 
fluids, although the ratio of the shear viscosity to entropy density is found to saturate the universal 
bound. Some of the retarded correlators are also computed in the high frequency limit for case study. 
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1 Introduction 

Gauge/gravity duality is a mighty device to study strongly interacting large N conformal field the- 
ories(CFTs) pQ-[5]. In this context, d-dimensional CFTs are mapped to gravity theories defined 
in asymptotically d+1 dimensional AdS spacetime, arising as its boundary theories. Especially, the 
fluid/gravity duality [6] -[9] has shed light on conformal fluid dynamics, which is low frequency (momentum) 
description of CFTs, and has provided much useful information such as transport coefficients of 
strongly coupled quark/gluon plasmas at RHIC and at the LHC [10]. Another interesting area is 
AdS/CMT [H]-[l5], which has recently witnessed a plethora of activities capturing interesting prop- 
erties of various condensed matter systems as thermal phase transitions between super conductors and 
normal conductors and their critical behaviors. 

While the initial interest in tangible examples in the subject concentrated on applications to AdS 
gravity theories and their CFT duals, the class of metrics of interest in gauge/gravity duality has 
been considerably enlarged including gravity systems, whose asymptotic behaviors are not AdS space 
time 0. The dual field theories defined on the asymptotically AdS^+i space is conformal and enjoys 
50(2, d— 1) conformal symmetry. However, it is widely accepted that one can generalize gauge/gravity 
duality prescription to the non-conformal field theory. Recently there is a deluge of research activities 
claiming that there are exemplar of dual field theories with gravity duals that do not possess conformal 
symmetries but another type of scaling symmetries, opulently found in Lifshitz geometry. 

For example, there are Lifshitz geometry, non-relativistic scaling metrics characterizing the so-called Schrodinger 
spacetime and even anisotropic space as dual geometries [3]- [47] in the literature. 
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From the gravity point of view, such violation of the hyperscaling symmetry [32]- [Hj can be 
represented in terms of the metric whose proper distance is not invariant under an appropriate scaling 
of the coordinates. It was shown in |38^l39j that the theory dual to the hyperscaling violating geometry 
can have a free energy including a logarithmic term in some parameter range, which takes a typical 
form leading to logarithmic violation of the area law in the presence of a Fermi surface. 

Especially, in a recent paper |48j . the authors have considered an interesting example of non- 
conformal fluid dynamics whose gravity dual is 4-dimensional Einstein-dilaton system with Liouville 
type potential V(4>) for the dilaton, parameterized by a constant rj as V ~ e^, where <fi is the dilaton 
field. It turns out that the dilaton solution has the form of a logarithmic function in r, where r is the 
radial coordinate of the bulk space time. This solution is asymptotically AdS when r/ = 0. However, it 
shows non-trivial asymptotics of the bulk metric as r — > oo when i] is non-zero, presenting a deviation 
from an asymptotically AdS space in terms of r/. Even if it shows such a non-trivial behavior, the 
boundary space time and thus its dual field theory on that are still well-defined. Moreover, one 
might expect that it may pose many sapid properties which are quite different from those in the usual 
conformal fluids and hence can be considered as another variant of fluid/gravity avatars(non-conformal 
fluid dynamics). 

In this note, we have studied some thermodynamic properties of the non-conformal fluid dynamics 
and computed transport coefficients to understand its properties more precisely. We have studied 
Hawking-Page phase transition of the system by using hard-wall model. It turns out that there is 
Hawking-Page phase transition at the critical temperature 

12 — 77 2 ±V ^=4 

which has a non-trivial ry dependence in it. ro is the location of the hard-wall and plays the role of an 
infrared cut off in the dual field theories. 

We have also resolved an issue concerning the curvature singularity in this background. The non- 
black brane geometry (black brane solution with zero size horizon, = 0, where is the horizon of 
the black brane) presents a curvature singularity in the interior, at r = 0, where i] ^ 0. However, in 
the hard wall model, it turns out that the singularity will be shielded by the hard wall and cannot 
be seen by an asymptotic observer. In those cases when the hard wall disappears, the black brane 
solution is always thermodynamically favorable than the non-black brane counter part. Therefore, the 
singularity will be located inside the black brane horizon. In summary, the curvature singularity will 
not be seen at all whether there is a hard- wall or not. 

We also investigated higher order corrections in u to the transport coefficients in hydrodynamic 
limit, where oo is frequency of the hydrodynamic modes in dual fluid system. For the shear viscosity 
described by the bulk gravitational perturbation h xy , we find that there is no emendation from 0(u) 
and the result thus obtained is consistent with the membrane paradigm result [59], where x,y are 
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boundary spatial coordinates. In addition to this, we have obtained the leading order corrections 
in small oj to the DC conductivity found in [48] governed by the transverse mode A y , where A y is 
y-component of the bulk U(l) fields. 

We have also computed the imaginary part of the Green's function using WKB approximation in 
the high momenta limit and found decay behavior of the imaginary part of the two point correlation 
function in the momentum space. We computed the diffusion constant from the bulk gravitational 
shear modes. The diffusion constant has a form of 



while the analogous counterpart for conformal fluids defined on 3-dimensional boundary is T> con f orma i = 
1/3. We also obtained the retarded Green's function from the gravitational shear modes and it shows 
poles at lo = —iDq 2 where q is the momentum of the hydrodynamic modes. 

The outline of the paper is as follows: In section [21 we study the Hawking-Page phase transition 
employing hard wall model and suggest a resolution of the curvature singularity problem in Einstein- 
dilaton theory. Section [3] is devoted to the computation of the leading order corrections in uj to the 
transport coefficients namely shear viscosity in section 13.1.11 and AC conductivity in section 13.1.21 
including a check of the universal bound on the ratio of shear viscosity to entropy density in the 
zero momentum limit. In section 13. 2\ we find a decay behavior of the imaginary part of the Green's 
function with momenta in high frequency and high momenta limit. Then in the following section I3~3l 
we present the crux of the gravitational perturbation with odd parity in the most general case and 
find the shear pole and the diffusion constant. Finally we conclude our work with some remarks in 
section [H 

2 Hawking-Page transition: A digression 
2.1 The Holographic setup 

We start with Einstein-Hilbert action with massless scalar field governed by Liouville type potential 
as 



where V(<j>) = 2Ae w \ A and T) are intrinsic parameters of the theory, G4 is 4-dimensional Planck 
constant and R is the curvature scalar. For convenience, we set G4 = 1. The equations of motion are 
given by 




(3) 




(4) 



(5) 
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The ansatz for the metric and dilaton solution are given by 

dr 2 

ds 2 = -a 2 (r)f(r)dP + + b\r) (dx 2 + dy 2 ) , (6) 

= 0(r). (7) 
Exact solutions of the equations of motion are obtained as 

o(r) = ^^iSe^V 1 , b(r) = b r ai , 
2^/12- r, 2 

/M-l-f^. (8) 



2?7 

and 0(r) = <fo - -j- — o ln ( r )' ( 9 ) 

4 + r] z 



where &o an d </>o are arbitrary integration constants, is horizon of the black brane and a\ and c are 
given by 



4 12 - 77 2 

A 9 aIld C = "1 ; 

4 + rf 4 + 7? 

For later convenience, we take bo = 1, cfio = and 



4(12 - r/ 2 ) 
(4 + r/ 2 



A = ~ A . JJ - (11) 



Such a choice of the parameters force ao = 1 and the Poincare symmetry of the boundary space time 
becomes manifest. The final form of the metric and the massless scalar are given by 

dr 2 

ds 2 = r 2ai [1 - (r/r fe H dr 2 + " + r 2 ^ (dx 2 + dy 2 ) , (12) 

r ZQl [1 — (r/r^J J 

= --^-ln(r), (13) 

where r is the Euclidean time obtained by wick rotation from the Lorentzian time t. We utilize this 
coordinate to argue thermodynamic properties of the Einstein-dilaton system, especially to discuss 
Hawking-Page transition in the next subsection. 

2.2 Hawking-Page phase transition 

In this section, we would like to discuss Hawking-Page transition of the Einstein-dilaton system. In 
some recent papers (e.g. [571 158]). the authors argued that there is Hawking-Page type of thermo- 
dynamic phase transition in the dual gauge theories by finding gravitational descriptions (identifying 
bulk on-shell actions with free energies in the dual field theories) where the dual gravity space time is 
capped off in the infrared near r = 0. This is the so-called hard-wall model. 

Notice that the radial cut off near r = in the dual gravity systems may correspond to the IR cut 
off in dual gauge theories. The theory dual to the non-AdS black brane describes a (2+l)-dimensional 
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relativistic non-conformal theory, which would come in handy as a toy model to understand the 
condensed matter system. If such system has an appropriate IR cut off (the size of the system), the 
dual gravity theory should also have an analogous counter part, a radial cut off. Such a cut off can 
be easily realized by introducing a hard-wall. 

In the following discussion, we will restrict ourselves in the regime as < rj 2 < 4, since it was 
already found in [38] that if rf > 4, then the black brane solution becomes unstable. To realize the 
Hawking-Page transition in this framework, we consider Euclidean version of the gravitational action 
of the Einstein-dilaton system ([3]) with an infrared cut off located at r = tq as 

i r r uv r/3 

Se = -— J dr J drd 2 x^ (R - 2(8<p) 2 - V(<j>)) , (14) 

where r\jy is the UV-cut off of the radial variable r to regularize the bulk gravity action. This is due 
to the fact that for both black brane and non-black brane backgrounds (the bulk geometry (|12p with 
r h = 0), Se is infinite, f is chosen as f = tq for non-black brane background since the spacetime will 
be capped off at that location. For the black brane background, f = Max(ro,rh), where l Max(A, B) 1 
means it will pick up the bigger value between A and B. j3 is the periodicity of Euclidean time, which 
will be specified in a moment. 

The precise integration can be performed for each of the cases. For non-black brane solution, the 
regularized bulk action is given by 

PV 3 



onbb _ P ^3 / 4+^2 4+^2 \ , . 

bE ~ 2^(4 + ^) yuv ~ r o I. U5) 



whereas the same for a black brane solution is obtained to be 



Se = 27r(4+V) ~ Max ^,r h )~^ j , (16) 

where /3b is periodicity of the Euclidean time for black brane geometry, which is chosen in such a way 
that there is no conical singularity at the black brane horizon. However, (3 is completely arbitrary. 
We compare the two Euclidean actions at the UV cut off, r = rjjv, where the local periodicity in the 
Euclidean time is the same, namely 



/3 



( r h \ 4+„2 

/3 B - (17) 



After all this manipulation, we compute the difference between the two bulk actions AS^, which 
is given by 

(18) 



ASe = km 

rjjy— >oo 



obb onbb 
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The precise computation shows that 



AS E = ^^r^ 2 for r >r h , (19) 

' 12 ~*>? 12->7 2 \ 

2r 4+ " - r^ +v for r < r h . 



47r(4 + r] 2 ) 

If ASV; is positive, then the non-black brane solution is thermodynamically preferable, otherwise black 
brane solution is preferred. Therefore, there will be a Hawking-Page type of thermodynamic phase 
transition at the critical temperature 

12 - ?7 2 4-^ 2 f=4 

±H ~ 4^(4 + 7/2) ' 1 j 

where we have used the expression of Hawking temperature of the Einstein dilaton system [48j . 



(12 -V 2 ) 
47r(4 + 77^ r ' 1 
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^ = i-rrr^r^ ■ (21) 



This result shows that for a fixed 77, the critical temperature increases when tq increases. More 
interesting is the behavior of the critical temperature Tfj which shows non-trivial r\ dependence for 
given tq. When r$ < ~ 0.0677, the critical temperature Tfj is a monotonically increasing function 
as i) runs from to 2, whereas when tq > 2l/9 1 e 2/ 3 ~ 0.4754, is monotonically decreasing. When 
< r o < 21/9^2/3 ; Tfj increases upto a turning point and after passing that point, it starts to 

(4+r, 2 ) 2 2(4+^ 2 ) 

decrease. This turning point is given by ro = 2 (>7 2 -i2) 2 e ^-12 in ro — plane. Some of the intuitive 
examples for these curves are provided in Figure. 1. As can be seen from Figure. 1, when 77 = 2, all the 
curves are converging to a single value, where critical temperature does not depend on the IR cut-off. 
In fact, this point is rather singular, since the black brane specific heat becomes infinite and Hawking 
temperature does not depend on the black brane horizon at all. 

Before closing this section, we note that "non-black brane geometry" contains a curvature singu- 
larity at r = 0. The curvature singularity in the interior of the "non-black brane geometry" is given 
by 

24(8 -r/ 2 ) 1 

(4 + r/ 2 ) 2 r 2^/(4+^)' [ ' 

which clearly shows a divergence unless rj = 0. 

However, this curvature singularity doesn't pose any problem to our previous discussion regarding 
the Hawking-Page transition since we have ripped off the singularity from the bulk space time by 
putting an IR cut off at r = ro. Even in the case of an infinitely small IR cut off ro, (as seen in the 
second line of (|19p ). the black brane solution is always thermodynamically favorable, therefore the 
singularity will be concealed by black brane horizon. Moreover, it turns out that the bulk on-shell 
Lagrangian density (free energy density in the dual field theory) is finite even for a zero size IR cut 
off, ro = when < ?? 2 < 4. In summary, the curvature singularity is not at all problematic. 
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H — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 
0.5 1 1.5 2 

Figure 1: This graph shows the behavior of the critical temperature, Tjj with different ro values. When 
ro = 1.10, Tjj is monotonically decreasing, whereas for ro = 0.05, Tfj is monotonically increasing. If 
ro = 0.30, Tfj increases upto r] ~ 1.1889, and then starts to decrease after passing that point. All the 
curves converge to a single point at t] = 2, where Tfj = 1/Att. 



3 Transport coefficients 

In this section, we will compute the shear viscosity, AC conductivity and diffusion constant %. 

9»» -> 9$ + V > A » 4 0) + SA n ■ ( 23 ) 

where gffi and are background metric and gauge fields, gffi is the background metric given by 
(|25p whereas is taken to be zero since there is no background gauge fields. For the later discussion, 
we will drop 5 from 5A^, then A^ will denote the perturbative bulk £7(1) gauge fields. 

In our calculations, we would like to use a new radial coordinate as u = — together with the 

' h 

rescaled variables t, x and y as {x,y} — > ^ {x, y}, t — > ^-t and choose 0o = ttt log(r/j)(the 
coordinates with tilde are primitive variables in metric([6])). Furthermore, we also take 



4(12 -r? 2 ^ 
(4 + r/ 2 ) 5 



A = ~ h I 2 , ( 24 ) 



as before [50]. Then, the background metric and dilaton solutions become 

du 2 

ds 2 = -a(u) 2 f(u)dt 2 + —- T — + a(u) 2 (dx 2 + dy 2 ) , (25) 

a{u) z f(u) 

27? 

■log(u), (26) 



where 



4 + r] 2 

a{u)=u a \ f{u) = l-u~ c . (27) 



2 We note that there are some of overlaps between our results in this section and those in [59]. Especially, for the shear 
viscosity and shear pole, computations are performed in the same background geometry and it shows quite consistent 
results with them. 
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Table 1: Modes Classifications of time dependent perturbations. 



Gauge Fields Each Decoupled Mode Global SO(2) Parity in x Parity in y 



NIL 


hxy 


Tensor 


Odd 


Odd 


NIL 


hxx hyy 


Tensor 


Even 


Even 


A x 


htx 


Vector 


Odd 


Even 


Ay 


h t y 


Vector 


Even 


Odd 


A t 


htt , hxx ~t~ hyy 


and 4> Scalar 


Even 


Even 


Table 2: Modes Classifications 


of frequency and momenta dependent perturbations. 




Gauge Fields 


Each Decoupled Mode 


Parity in y 






Ay 


h ty and h xy 


Odd 






At and A x 


htt, h xx , hyy, htx and 


Even 





In the following discussions of the transport coefficients, we will mostly use the rescaled coordinate. 
In this coordinate, the horizon is located at u = 1. The black brane temperature is given by 

T„ = ±, (28) 

and entropy of the black brane is S = where V2 = J dxdy is the spatial volume of the boundary 
space time. 

To obtain hydrodynamic modes in the dual fluids, we turn on metric perturbations together with 
dilaton perturbations. For a better computation, we classify the gravitational and dilaton modes as 
follows. The background metric solution ([25]) enjoys 5*0(2) global symmetry (spatial rotation in x — y 
plane) together with parity symmetries as x — > —x and y — > —y. Therefore, perturbations can be 
classified in accord with these global symmetries when they depend on time only. They are given by 
Table. 1 

However, once one turns on spatial momenta, it breaks these symmetries into smaller symmetry 
group. Without any loss of generality, we can choose the direction of the spatial momenta along 
x-direction. In this case, global 50(2) symmetry breaks into U(l) and only y — > —y parity symmetry 
is retained. Global £7(1) is just a translational symmetry along y direction and is not relevant for the 
mode classification. Therefore, we sort out the perturbative modes by y parity symmetry. They are 
given by Table. 2. 

For the shear viscosity one considers a tensor mode of the metric fluctuation denoted by h xy . 
For the computation of electrical conductivity using Kubo's formula, it is sufficient to consider gauge 
fluctuations around the black brane solution and keep other fields unperturbed. In order to determine 
diffusion constant, it is enough to consider perturbations in (ty) and (xy) component of the metric 
tensor. 
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3.1 The shear viscosity and AC conductivity in the zero momentum limit 

In the zero momentum limit, the gauge fluctuations are decoupled from the metric ones in Table 2. In 
this limit, the study of the transverse gauge fluctuation with nonzero frequency provides information 
about the AC conductivity, which reduces to the DC conductivity in the zero frequency limit [56]. In 
addition, we can find the shear viscosity from the bulk gravitational perturbation h xy . 

3.1.1 The shear viscosity 

Consider a tensor mode of the metric fluctuations, h xy . In the zero momentum limit, using the 
following Fourier transform as 



tfM = / s-e-^u,"), (29) 



dt 

Einstein equations for h xy at the linear order in weak field expansion is given by 



' tl An \ up' 

o = K" + [j + v) < + ^72^' (30) 

where the prime denotes the derivative with respect to u and hy = hy(u, uj). The solution of the above 
equation can be expanded as 

h x y = f {G (u) + wGi(u) + uj 2 G 2 (u) ■■■), (31) 

where 5 is a constant. Notice that the horizon lies at u = 1 and / = at the horizon. Due to the 
singular structure of (|30p at the horizon, the near horizon behavior of hy is given by 

h x y = f ±lu > /c , (32) 

where f- luJ / c (f+ luJ / c ) implies the incoming (outgoing) solution. In order to study the transport coef- 
ficients from boundary retarded Green's functions, we should choose an incoming solution. We solve 
the equation (|30p order by order in small frequency and the zeroth order equation in uj is 

0=(« 4si /GoM')'- (33) 

The regular solution(at the horizon) of the equation is given by Gq(u) = cq where cq is an arbitrary 
constant. The first order equation in uj is of the same form as (|33p and thus its solution is given by 

16 log u - (4 + if) log (u - u 16 /( 4+? ? 2 ) 



Gi(«) = ci + c 2 £ '-. (34) 

12 — r/ z 

The regularity at the black brane horizon forces c 2 = 0. On AdS boundary, we impose (the so-called) 
"vanishing condition" as 

h x y (u -> 1) = c , (35) 
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then we set c\ = under such boundary condition. 

It may be interesting to ask how the form of the shear viscosity gets modified after turning on 
higher order corrections in small frequency. Such corrections can be obtained if we take a parametric 
regime where frequency is small but finite. Since the first order correction in u vanishes, G\(u) = 0, 
in order to study non-trivial small frequency corrections, we solve the second order equation in ui in 
((301) . The second order equation in u) provides a solution G2(u), which turned out to be of an integral 
form as 



G 2 (u) 



dv! 



w 



/4ai 




„Co(l-U *+v r ) 



It" 



II— ^ 

u 4 +1 



+ c 3 + c 4 , 



(36) 



where C3 and C4 are integration constants. C3 is determined from regularity condition at the horizon and 
it turns out that C3 = 0. The vanishing condition at the black brane horizon fl35f) forces G*2(u — > 1) = 0. 
This condition gives C4 = 0. In fact, analytic form of G' 2 (u) can be obtained as 



G' 2 (u) = 



co 



4ai 



/(«) 



4 + rf 



12 



V 



;i 



4 + if 



12 



12-17 

4+r, 2 



4 + 7?^ 



12 



■Zn(l 



12-?7_ 
4+I7 2 



(37) 

where is Hyper geometric function. 

Using all the above information, we evaluate bulk action upto equations of motion for hy, which 
is given by 



S 



1 

'32?r 
1 

'32tt 



d 3 x u^f(u)h x y (u) h*'(u) + Sct 



(38) 



dujd 2 qcQ bJ Co. 



ICO — CO 



4 + 7T 

12-77 



-7 



PG 



4 + r) 2 
12 - if 



+ O(co 3 ) 



S c t is counter term action 0. As long as < rj 2 < 4, PG I 



where cq is the boundary ^value of /i^(at u = 00), 7 is Euler number, PG is Polygamma function and 

is purely real. The retarded Green's 

4 + ?? 2 



12-7?^ 



function can be read off from (1381) as 
G 



r xy,xy 



1 



id 



4 + T7 2 



12 



PG 



12 



which shows corrections in 0{ui 2 ). We define shear viscosity function as 



xy,xy \ 



(40) 



(41) 



and in the uj = limit it is identical to the shear viscosity obtained from the Kubo formula. We 
evaluate r)(oj) in the hydrodynamic limit and it may show sub-leading corrections in uo to the shear 



3 Under the vanishing condition, the boundary value of /i^(at u = 00) is not precisely Co, but get some corrections. 
4 The boundary value of hy is given by 



co 



dv! 



u' 4 ^f{u') 



II. 



(ill 



„c (l - U ) 



(39) 
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viscosity obtained from the Kubo formula. However, it turns out that it does not contain higher 
order corrections in small frequency upto 0(ui), since the shear viscosity function cares only about the 
imaginary parts of the retarded Green's function. Therefore, 

= tL + ° {uj2) - (42) 

The shear viscosity saturates the universal bound upto the first sub-leading order correction in u, 
s = lb> being consistent with the result obtained from the membrane paradigm [49} 150] . where s is 
the entropy density in the rescaled coordinate, s = |. In addition, our study on the higher order 
corrections to the shear viscosity shows = in the zero frequency limit. 

3.1.2 AC conductivity 

In [481150] . DC conductivities of the non-conformal holographic fluids was investigated. In this section, 
we will study AC conductivity of the fluids by using its gravity dual, Einstein-dilaton system. 
The gauge fluctuation coupled with dilaton field is given by 

SA = ~I^j ^V^ge^FTu*, (43) 

where g 2 is the gauge coupling. The equations of motion of transverse modes, A x and A y in the zero 
momentum limit (See Table. 1), is given by 



where A indicates either A x or A y and 



b= 2 ±^>. (45) 
4 + r] 2 



Notice that in the zero momentum limit, two transverse modes A x and A y satisfy the same differential 
equation due to the global SO(2) symmetry(rotation in x — y plane). 
The near horizon solution of A has a form 

A ~ f ±luJ ' c , (46) 

in which the positive (negative) sign is for the outgoing (incoming) boundary condition at the horizon. 
Imposing the incoming boundary condition on A, the solution in the hydrodynamic regime (co « 1) 
can be expanded in small frequency as 

A = r iul,c (Go(u) + wGi(u) + oj 2 G 2 (u) + •••), (47) 

in which the ellipsis indicates the higher order corrections in co. 
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The zeroth order equation in uj is given by 

= ufG , > + (bf + uf')G> , 

and its solution is obtained as 



G = ci + 



° 2 „l-&+c p 
It 2-^1 



l-b + c 



l-b + c l-b+c ( 
,1,H ,« 



The regularity on Go at the horizon forces C2 to be zero, then we get 

G = c\. 

The first order equation in u is 



1i(li c — 1) 



This equation has the following solution: 



Gi = c 3 + 



C4 oA-b+c P 
U 2-T1 



1 - 6 + c 



1-6 + c 

At the horizon, the regularity fixes C4 to be 



,1,1 + 



u 2 (u c - 1) 
1-6 + c 



ic\ u c - 1 
w H log — 



C4 = ZCi 

and the vanishing condition of Gi determines C3 as 

ici [£G + iir + PG (0, ±=*±£)1 
c 3 = L - i 

c 

where EG and PG imply the Eulergamma and Polygamma functions respectively. 
Finally, the second order equation in uj is obtained as 

= 8 U (u b fG' 2 ) + £L u -2-4a 1+ fe-2c (n 2+2 C _ ^ j + ^4a 1+ c [(1 _ fe + ^ _ 3^] . 

Its solution has an integral form like 

G 2 = - 



" /1 + /2 + C5 , 
au — h c 6 , 



where 
/1 = 



ci 



ci 



2jFi (_ l=|t£,l,- 1=6, ««) u l-4ai+6+c ( l-4al+iH-c >1) l-4al+6+2c ^c) 



(l-b + c) u 1 -^ 



+ 



(3 + 8 ai - b + 3c) n 1+4al+c 2 -Fi (l, 



1 - 4al + 6 + c 



l+4al+c o 1 l+4al „,c^ 



;,2 + i±±ai,u c ) 



(1 + 4oi + c) 2 

9 u 2+4ai-b+2c p / 2+4ai-fr+2c i 2+4ai~fe+3c c 
2 11 c , , c , 



2 + 4ai - 6 + 2c 



+ 



+ 



i 2 + 4 °i- 6 + 2c 2^(1,1^,2 + 1^,^) (1-6 + c) / 2 + 4 ai 

*0 [//.__ ~\~ 



1 + 4ai + c 
(l-6 + c)u 1+4ai+c 

; r < ITT 

c (1 + 4ai + c) 



+ ?£Af 



c (1 + 4ai + c) 

l-6\ u c -l 
+ log — 
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HJ\f and B denote the harmonic number and beta function respectively. C5 and c% are integration 
constants, which can be determined in terms of c\. To determine them, after performing the integration 
in (|56j) near the horizon, we impose the regularity and vanishing condition. The regularity at the 
horizon fixes C5 to be 



c (1 + 4ai + c) [ V c / V c 

l-4ai+6 + c\ / -1 + 6 

+-{ 3, ) -PG[ 0,-1 + 



(58) 



where the vanishing condition sets cq = because the integration in (|56|) does not give any constant 
at the horizon. It is worth noting that the fact cq = implies that the oj 2 order solution does not 
contribute to the boundary value of A. 

Imposing the Dirichlet boundary condition at the asymptotic boundary like 

A = lim A(u). (59) 

u— >oo 

c\ can be determined by 

c cAq , 6Q s 

Cl c-TTco + iuj [EG + PG (0,1=2*)]' 

Bulk action (|43]) together with (|59p and bulk equations of motion becomes 

S A = - lim I d 3 x A u b f(u) A(uY. (61) 

u— >oc q y 

One can see that the terms proportional to u~ b in A{u)' can only contribute to the finite part of the 
retarded Green's function of the transverse mode. Using this fact, the imaginary part of the retarded 
Green's function reads as 

G m = %(l + ?- U ). (62) 



g 2 \ c 



After all dusts get settled, the AC conductivity becomes 

_ G xx 1 / (4 + r? 2 )^ \ 
(J AC = — = ~z 1 + — = — w . (63) 

In the zero frequency limit (co —> 0), the AC conductivity reduces to the DC conductivity, which is 
consistent with the result obtained in [38] and also in [50j by the membrane paradigm. By using 
u = r/rh and u) = r/jW, we get 

aAc = ?i 12-, 2 J r ^ + 4^{-i2^r) h (64) 

where the mass dimension of u is the same as that of In the hydrodynamic limit, the conductivity 
increases with increasing frequency. The conductivity increases (or decreases) with temperature for 
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a < (or a > 0), which shows that for a > the dual field theory has a metal- like behavior 
(DC conductivity decreases with increasing temperature), whereas for a < it shows electrolyte- 
like behavior with a decreasing DC conductivity with temperature. Especially, the known metallic 



behavior for the resistivity occurs at a 
to the temperature p ~ Tjj. 



27) 



which shows the resistivity to be linearly proportional 



3.2 Green's function at high momenta 

In this subsection, we will consider the high momenta (lj, \q\ >> 1) domain of gauge perturbations 
discussed in the previous subsections. We use a similar technique as given in |51j . We first turn on 
the momentum q along x direction. Then the gauge fluctuations At and A x are coupled while A y 
represents the transverse mode. This transverse mode satisfies 



= A» + [ b - + QA' + ^lA 



u f 



u 



4aij2 



(65) 



Now, we make the transformation 



Ay(u) 



b \ 1/2 



U c - 1 



Ay(u). 



to bring (J65J) in the Schrodinger form 



A'^u)= [uj 2 Q{u) + R{u)] A y {u), 



where 



and 



<?(») 



KM = ^ 



+ 



P 2 f-1 

M 4aij2 

2c(c-6 + l) 



b{b - 2) 



(66) 

(67) 

(68) 
(69) 



(u c - l) 2 ' u c -l 
with p = \q\/uj. 

In the large frequency limit (lj » 1), the term proportional to lj 2 dominate the potential and the 
WKB approximation is possible. To apply the WKB analysis, we first expand A y (u) as 



A y (u) = exp 



i=0 



(70) 



After substituting this WKB ansatz in (f6T|) we get, at zeroth order in 1/lj 

So = ± J VW) 



(71) 



For the spacelike momenta, lj 2 < \q\ 2 (p > 1), the potential Q(u) is positive in the range uq = 
(1 - \)~ 1,c < u < oo while it is negative otherwise. Hence, the WKB solution for (|67p decays 
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exponentially in the interval {uq,oo) and oscillates in the interval (1,uq). Physically, this implies 
that the incoming wave first tunnel from u = oo to u = uq before it reaches the horizon u = 1. 
The imaginary part of the retarded Green's function G yy , which is proportional to the tunneling 
probability, is then given by 



SsG yy (u,q) ~ exp 
and in the limit where to <C \q\, it becomes 



(72) 



QGy y (uj,q) ~ exp 



where 



A 



-2up 



u 



-2ai 



1 y/l — U c . 

20Fr[(2 Ol - l)/c 



exp 



A|< 



(73) 



(74) 



_cr[(4ai -2 + c)/2c] 

Note that, the above expression for Green's function is independent of dilaton parameter a. In other 
words, in the regime of high momenta the dilaton field does not affect the retarded Green's function. 
This decaying nature of the imaginary part of the Green's function is similar to the one obtained in 
[51] . However, in our case the decaying probability depends upon the parameter 77. The constant A is 
a positive number, bounded from below by its minimum value as A > \ m in = "irrg^sT^ ~ 2.8044 in 
the region < rf < 4, which ensures a decay behavior of the Green's function in such window for r]. 

3.3 Gravitational perturbations with odd parity in y and Retarded Green's func- 
tion 

In this section, we will turn on the parity odd modes, ht y and h xy and obtain shear modes in the dual 
fluid dynamics. Non trivial equations of motion up to first order in weak field expansion are given by 



W, 



Wi 



ty 



W, 



uX' 1 (v)+qf(v)X' 2 (u)=0, 

X l W + i \ X lW - 4 , w7l - 4/ rj7 \ = °> 
a(u) a\u)f(u) a^(u)f{u) 

a'{u) i f'(u)\ t u 2 X 2 (u) quX^u) 



A 2 (u) + X 2 (u) 4 - + — + - 
V a(u) f(u)J a 



0. 



(75) 
(76) 

(77) 



x*(u)p(u) a\u)p{u) 

where the perturbative fields are obtained in the frequency(momentum) space by the Fourier transform 
defined in the previous section and the fields X\(u) and X 2 (u) are given by 



h" y q = a 2 {u)Xi(u) and h^ y q = a 2 {u)X 2 . 



(78) 



X\ and X 2 definitely have frequency and momentum dependence but we suppress the explicit expres- 
sion of those for the sake of calculation. To solve the above set of coupled differential equations, we 
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rewrite them only in terms of X\ by combining them (with the aid of a little algebra) [511 [52] and the 
final equation to solve is given by 

= (/(«) (a\u)X[{u))') - (V - X[(u). (79) 

This is a third order equation and we try the following ansatz for the solution: 



( 12-«7 \ (-A)(4+r, 2 ) 
l- u G(u), (80) 

where 

oo 

G{u) = wVXii(u), (81) 

C is an overall constant and the factor in front of G(u) is for the ingoing boundary condition near 
black brane horizon. It turns out that each of Y\ii are given by 

*i,oo = u'^ 1 (82) 

/ 12-») 2 

16 4i(r? 2 - 12) is /l-u 

n,io = C> 3+^ + l' tA ■ ^[ 2 u ^ I du\ — r | (83) 



-A(4 + r? 2 ) 2 



1 — u 4+ri 



— 



° loU + -A(4 + r? 2 ) 2 12 _ r?2 >i 12-^2^ 

4 + r? 2 _i2-^ 2 
12 - r? 2 v 7 

4u 
-A(4 + r? 2 



— / 1 4fcu, \ 

Yifii = 0, Fi )02 = u { Q 2 + A ^ ^ 2 . ) , and so on, (84) 



where Cfj are O(l) constants and 2-H1 is Hyper geometric function. To obtain shear pole (diffusion 
constant), we plug the solutions into (|76|) . With boundary values of the fields, X\(u) and ^(it) as 

X^u-^oo) =X X and X 2 (u -»• oo) = X 2 , (85) 

we obtained the overall coefficient C in terms of those as 



-A(4 + ?] 2 ) / g 2 Xi + gwX 2 

4+»? 2 

From this expression, we ensure that the pole of the green's functions from such modes (odd parity 
modes) will be given by 

u = -W<? and V = i±4 = (87) 
y 12 — 77 2 1671-7^(12 - r] 2 ) v ; 

where T> is the diffusion constant. Diffusion constant is usually given by T> = fj/(e + P), where fj is 
the shear viscosity, e is the energy density and P is the pressure of the dual fluid system. For the 
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fluids whose gravity dual is long wavelength excitation in AdSd+i spacetime, T> = 47r y |53j (See the 
table. 1 in this reference). By using the definition of the cosmological constant A in (|24p . one realizes 



that T> = inT for our case too. 



In fact, the relevant parts of the bulk on-shell action to obtain the retarded Green's functions in 
this case is given by 



'bulk— os 



r 2 h ai f-A(A + rj' 2 
32^ 



dujd 2 q 



and from this expression, the retarded Green's functions can be read off as 

8 4+17 2 

R (47rTff)^ ( A + rf 

Gr 



16vr ~\\2-tP) q 2 -iu^ 

8 4+17 2 

R (AirT H )^ / 4 + r? 2 \ *=V u 2 

xy ' xy 167T \12- V 2 J [ j 



4+r/ 2 



i+v 2 



4 Discussion 

In this note, we have discussed at length the transport coefficients of the non-conformal fluid system, 
whose gravity dual is Einstein-dilaton theory with a Liouville type dilaton potential parameterized by 
an intrinsic constant r]. The ratio of the shear viscosity to entropy density in zero frequency limit turns 
out to saturate the universal bound. However, the other transport coefficients in the hydrodynamic 
limit such as AC conductivity of the i?-charge current, diffusion constant of the shear modes of the 
non-conformal fluids depend on r\ non-trivially. 

We also computed the retarded Green's functions of R-charge current in the high frequency limit. 
It turns out that the retarded Green's functions show exponentially decaying behaviors in momentum 
space as G R ~ e~ x ^ q ^ T , where the constant A depends on tj in a non-trivial way. 

For our next project, we will consider higher order conductivities from bulk £7(1) gauge fields in 
both frequency and spatial momenta and examine their dependence on ??. We will also try to apply 
the WKB approximation scheme to other bulk perturbations in the high momenta limit and check its 
decaying behaviors. So far the interpretation of the retarded Green's functions in the high momenta 
limit is not clear in the dual fluid dynamics. We would like to shed some light on this high momenta 
regime by studying the Einstein-dilaton system in gory details. 
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